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Abstract 

Neutron scattering was used to determine the spin structure and the magnon 
spectrum of the Mott-Hubbard insulator YTiOs. The magnetic structure is 
complex, comprising substantial G-type and A-type antiferromagnetic compo- 
nents in addition to the predominant ferromagnetic component. The magnon 
spectrum, on the other hand, is gapless and nearly isotropic. We show that 
these findings are inconsistent with the orbitally ordered states thus far pro- 
posed for YTiOa and discuss general implications for a theoretical description 
of exchange interactions in orbitally degenerate systems. 
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The interplay between spin and orbital degrees of freedom in transition metal oxides 
has been a subject of investigation since the 1950's. As one of the salient outcomes of this 
effort, the "Goodenough-Kanamori rules" provide a description of the exchange interactions 
between magnetic atoms, and hence the magnetic ordering pattern, in terms of the relative 
orientation of valence orbitals on neighboring lattice sites. This field has recently moved 
back into the center of attention as advances in materials preparation have made it 
possible to investigate not only the static spin and orbital arrangements, but also the spin 
and orbital dynamics in (i-electron systems. In cubic manganites, the focus of much of the 
recent attention, the spin wave excitations have been studied extensively 0, and excitations 
associated with the cZ-orbital degrees of freedom have recently also been detected [0]. 

The cubic crystal field in the perovskite structure splits the degenerate (i-orbital manifold 
of a free transition metal ion into a lower-lying triplet of t2g symmetry and a higher-lying Cg 
doublet. In the manganites, where the Cg doublet is partially occupied, coupling to the lattice 
through the Jahn- Teller effect lifts the orbital degeneracy and generally pushes the orbital 
excitations to energies much larger than the magnon band width. The spin and orbital dy- 
namics are thus largely decoupled, and quantum effects are suppressed. For t2g orbitals, on 
the other hand, the higher degeneracy and the more isotropic, less bond-directional charge 
distribution reduces the lattice coupling, and one may expect a more dramatic interplay 
between the orbital and the spin dynamics. The manifestations of this interplay should be 
most apparent for a single d-electron in the t2g manifold, a situation realized in the pseudocu- 
bic titanates. Indeed, the magnetic properties of the Mott-Hubbard insulator LaTiOa are 
difficult to reconcile with predictions based on the standard Goodenough-Kanamori rules. 
The G-type antiferromagnetic ground state ||^ of LaTiOs is in conflict with the predictions 
of all electronic structure calculations thus far reported Further, its ordered moment 

is smaller and its spin wave spectrum is more isotropic than predicted by conventional 
superexchange models 0. 

The theoretical implications of the unusual magnetic properties of LaTiOa have remained 
controversial. The theories thus far proposed invoke disparate effects ranging from a subtle, 
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hitherto unobserved lattice distortion to a novel spin-orbital resonance PJTO[]. In order 



to extend the empirical basis for these model calculations, we have investigated the micro- 
scopic magnetic properties of YTiOs, a sister compound of LaTiOa whose larger 0-Ti-O 
bond angle results in a reduced electronic band width and an increased Mott-Hubbard gap 
rT|,|T2|. In contrast to LaTiOa, the ferromagnetic ground state of YTiOa (Ref. is in 



accord with electronic structure calculations ||,^ . We have used neutron diffraction to show 
that the full magnetic structure of YTiOa, while predominantly ferromagnetic, is actually 
noncoUinear and hence more complex than previously assumed. Further, an inelastic neu- 
tron scattering study reveals a magnon spectrum that cannot be explained in terms of the 
orbital ordering patterns thus far proposed. We discuss our results in the light of current 
theories of superexchange in the presence of orbital degeneracy. 

The sample was an untwinned single crystal of volume 0.4 cm^ grown by the floating zone 



technique as described elsewhere [|TT|. The neutron scattering experiments were performed 
at the IN22 triple-axis spectrometer at the Institut Laue-Langevin in Grenoble, France, and 
at the E5 four-circle diffractometer at the BER-II reactor of the Hahn-Meitner Institute in 
Berlin, Germany. For the measurements on IN22, a pyrolytic graphite (PG) monochromator 
and a PG-analyzer, both horizontally collimating, were used. Depending on kinematic 
constraints and resolution requirements, the final wavevector was fixed at 1.64 A~^, 2.66 
or 4.10 A~^. On E5, neutron wavelengths 2.36 A or 0.902 A were selected by a PG or a Cu 
monochromator, respectively, and the data were collected with a two-dimensional position 
sensitive ^He-detector. 



A set of 1359 structural Bragg reflections taken on E5 [Tj| could be refined with the 
orthorhombic space group Pnma and room temperature lattice parameters a = 5.3584(9) A, 
b = 5.6956(8) A and c = 7.6371(11) A, in agreement with previous reports @|1^. In order 
to simplify the comparison to model calculations, we sometimes use the pseudocubic (c) 
unit cell with lattice parameters a' = a/^/2, b' = and c' = c/2 instead of the full 

orthorhombic (o) cell. Our data provide no indication of a structural phase transition below 
room temperature. 



The presence of magnetic intensity could be established from the temperature depen- 
dence of the primitive Bragg reflections at low momentum transfer, some of which show 
a spontaneous increase in intensity below Tc = 27 K (Fig. |T]). The reflections (2, 0, 0)o, 
(0, 2, 0)o, (1, 1, 0)o and (1, 1, 2)o show the strongest magnetic intensities, indicating a 
predominantly ferromagnetic coupling of the Ti-moments as previously reported |^. In 
addition, however, magnetic Bragg reflections were also observed at (0, 1, l)o, (1, 0, l)o, 
(0, 0, l)o and (0, 0, 3)o (Fig. |l]b). This indicates the presence of hitherto unobserved 
staggered components of the ordered moment. A full analysis of the magnetic diffraction 
pattern at T = 10 K gives the following components of the ordered moment per Ti'^+ ion: 
a ferromagnetic component of 0.544(11) fis along c; a G-type antiferromagnetic component 
of 0.082(12) fiB along a; and an A-type antiferromagnetic component of 0.047(13) fiB along 
b. In agreement with Bertaut's representation analysis [|16|, the magnetic structure can thus 



be described by the basis functions Gx, Ay and Fz |T^. A pictorial representation is given 
in the inset to Fig. 0. The total magnetic moment at 10 K is 0.553(11) yU^/Ti^^-ion. Since 
the magnetic Bragg peak intensity is not saturated at 10 K, this extrapolates to 0.715 fiB 
at T = 0, somewhat lower than the saturated moment of 0.84 fiB previously reported for 



samples with proportionally higher Curie temperatures |TB[ . 

We now turn to the determination of the spin wave spectrum by inelastic neutron scat- 
tering. Fig. 1^ shows representative constant-q scans at T = 5 K and 50 K, that is, above 
and below Tc- As shown in the figure, magnons can be distinguished from phonons on the 
basis of their temperature dependence. While the former become overdamped above Tc, the 
latter only show a slight intensity change due to the thermal population factor. Since the 
lineshape of the observed features is strongly influenced by the spectrometer resolution, a 
deconvolution is required to accurately extract the spin wave peak positions. We have used 



both the standard Cooper-Nathans procedure and a Monte Carlo ray-tracing routine p!9 
Some of the resulting profiles are shown in Fig. |^. 

The inelastic neutron scattering data are summarized in Fig. ^. Some features of the 
spectrum are immediately apparent. First, the spectrum is gap less: by high resolution 



constant-q scans at the zone center we have obtained an upper bound of 0.3 meV on the spin 
wave gap. Second, the spectrum is almost isotropic: a surprisingly good fit can be obtained 
by an isotropic Heisenberg model which gives a dispersion of the form huj = 65171(1 — 7q) 
with 7q = ^{cos{qxa') + cos{qyb') + cos{qzc')) and S = 1/2. The best fit of the spectrum 
is obtained by a nearest-neighbor superexchange parameter J = —3.0 meV, corresponding 
2^ to Tc= [1A5S{S + 1) - 0.18)]\J\/kB = 32 K, which agrees well with the experimental 



value. A small systematic deviation from the result of this fit is noticeable, albeit within 
the experimental error bars (Fig. 1). The fit can be improved by introducing either a 
subtle (~ 6%) difference between the exchange parameters in the a6-plane and along the 
c-axis, or by incorporating a small exchange anisotropy that is also required to explain the 
noncoUinear spin structure. Specifically, the leading anisotropy term in the spin Hamiltonian 
of the titanates is of cubic symmetry ||10| and modifies the magnon dispersion as follows: 



huj = \/3| J|(l - 7q) + A + A{1 - cos{q,a')) x v'3| J|(l - 7^) + A + A{1 - cos(g^6')) , (1) 

where A is the anisotropy parameter and A = 0.093A^/| J| is the spin wave gap. The 
best fit gives A = 0.8 meV and J = —2.75 meV so that A = 0.02 meV, consistent with 
the experimental upper bound. The cubic anisotropy also induces spin canting in a pattern 
identical to the tilting pattern of the TiOg octahedra in the Pnma structure. This is 
indeed the dominant component of the observed canting pattern (Fig. [l|). An isotropic 
Heisenberg Hamiltonian with a small exchange anisotropy thus provides a good description 
of the magnetic ground state and dynamics. 

While the neutron scattering data can be fitted by single, resolution limited profiles over 
most of the Brillouin zone (Fig. ^), the profiles near the (|, |, |)c point, where the magnon 
crosses an optical phonon branch at 19 meV, are broader than the instrumental resolution 
(Fig. 0b). The same effective broadening was also observed in the spin-flip channel of 
a polarized-beam experiment. The best overall fit of the constant-Q scans near the zone 
boundary is obtained by two-peak profiles corresponding to two sharp but nondegenerate 
magnons (Fig. 2b). The splitting of the magnon degeneracy probably originates from a 



hybridization of the magnon with the optical phonon. Magnon-phonon hybridization has 
been observed in other transition metal compounds and was attributed to a modulation of 
the crystal field potential by the phonon ||21]|. However, as the two peaks are not completely 
resolved, other origins (such as a finite magnon lifetime) cannot be ruled out entirely at 
present. 

Leaving these subtleties aside, we now discuss the overall features of the magnon spec- 
trum in terms of the orbitally ordered state predicted by band structure calculations 
and reported to be consistent with NMR and neutron [^] form factor data. In this 
state, the (i-electron occupies the following orbitals on the four inequivalent Ti sites of the 
Pnma structure 

|V')l,3 = Ci\dy^) ± C2|4y), 

1^)2,4 = Ci\d^^) ± C2|4s,)- (2) 

Using the superexchange Hamiltonian given in Ref. |^ and the above wave functions, 
the expectation values of the c-axis and in-plane exchange couplings are: 

Jc = JsElin + r2r3)(l - n^y) - (ri - r2)](l - n^y)/2, 
Jab = JsEiin + r2r3)nly - (ri - r2)(l + n^y) /2]/2. (3) 

Here, Jse = /U (where t stands for the Ti-Ti hopping parameter and the Kanamori 
parameter U for the Coulomb repulsion on the same t2g orbital) represents the overall 
superexchange energy scale, and n^y = is the fraction of the d^y orbital controlling the 
orbital state. The coefficients ri = 1/(1 — 3//), r2 = 1/(1 — r]) and = 1/(1 + 2//) originate 
from the Hund's rule splitting of the excited ^23 niultiplet hj rj = Jh/U. The parameters 
Uxy and rj control the exchange couplings along the different axes and their ratio, as shown 
in Fig. g. 

In principle, one may obtain isotropic exchange parameters (that is, Jc = Jab) of ferro- 
magnetic sign ||2^ by using an orbital state with n^^y ~ 0.6 (close to the value reported in 
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Ref. |2^). However, this requires a large, unrealistic value of the Hund coupling. For repre- 
sentative values of Jh ~ 0.64 eV and U — ~ 4 eV, [^j one can estimate 77 ~ 0.12 which 
is inconsistent with isotropic, ferromagnetic exchange parameters (Fig. Moreover, an 
unphysical fine tuning of the orbital state is required to obtain a spin Hamiltonian of cubic 
symmetry. Because of the high sensitivity of the spin interactions to the orbital state, even 
a small (~ 5%) change of n^y leads to a strong spatial anisotropy of the coupling constants 
and may even reverse their relative sign (inset in Fig. 

An explanation of the small spin wave gap A < 0.3 meV presents a further serious 
difficulty for the conventional Goodenough-Kanamori picture. The TiOg octahedra of YTiOs 
are elongated, leaving four almost equal Ti-0 bonds To a first approximation, the 

crystal-field ground state should therefore be a quadruplet with substantial unquenched 
orbital angular momentum and a large magnon gap. Using the wave functions of Eq. |^ with 
parameters that result in = Jab ~ — O.IJ^^;, we have obtained a quantitative estimate of 
A Jse{^/ DY the gap due to the symmetric exchange anisotropy, where A ~ 15 — 20 
meV is the spin-orbit coupling constant and D is the crystal-field splitting of the two lowest- 
lying Kramers doublets. In addition, the Dzyaloshinskii-Moriya interaction cants the spins 
by an angle ~ (J5£;/3| J|)(A/-D) away from the c-axis and further increases the gap by 
I J|0^. To be consistent with the experimentally determined canting angle and upper bound 
for the gap, D must be a least 200 meV which is hard to reconcile with the tetragonal crystal 
field environment of the Ti^"^ ion. 

These considerations imply the need to go beyond the standard Goodenough-Kanamori 
approach in order to arrive at a theoretical description of the nearly isotropic magnon spec- 
trum of YTiOa. While the experimental situation p2| , p3| precludes an orbitally disordered 



state such as the one proposed for LaTiOa, ||^|T^ quantum zero-point fiuctuations in the or- 
bital sector naturally lead to magnon spectra with diminished spatial and spin anisotropics 
even if orbital order is not entirely obliterated. This has been substantiated in a recent 
explicit model calculation . 



In summary, we have shown that the magnon dispersions of YTiOa are an extremely sen- 



sitive gauge of the orbital state. The measured dispersions are inconsistent with the orbitally 
ordered states thus far proposed and point to the importance of orbital zero-point fluctua- 
tions that go qualitatively beyond the conventional Goodenough-Kanamori picture. YTiOs 
is an excellent model system for a full quantum many body theory of the superexchange 
interaction in orbitally degenerate systems, an interesting subject of future research. 



8 



REFERENCES 

[1] For a recent review, see Y. Tokura and N. Nagaosa, Science 288, 462 (2000). 

[2] See, e.g., J.W. Lynn et a/., Phys. Rev. Lett. 76, 4046 (1996); F. Moussa et a/., Phys. Rev. 
B 54, 15149 (1996); K. Hirota et al, J. Phys. Soc. Jpn. 65, 3736 (1996). 

[3] M. Saitoh et al, Nature 410, 180 (2001). 

[4] D.A. MacLean, H.N. Ng, and J.E. Greedan, J. Sohd State Chem. 30, 35 (1979). 

[5] T. Mizokawa and A. Fujimori, Phys. Rev. B 51, 12880 (1995); ibid. 54, 5368 (1996); T. 
Mizokawa, D. I. Khomskii, and G. A. Sawatzky, ibid. 60, 7309 (1999). 

[6] H. Sawada and K. Terakura, Phys. Rev. B 58, 6831 (1998). 

[7] B. Keimer et al, Phys. Rev. Lett. 85, 3946 (2000). 

[8] M. Mochizuki and M. Imada, J. Phys. Soc. Jpn. 70, 2872 (2001). 

[9] G. KhahuUin and S. Maekawa, Phys. Rev. Lett. 85, 3950 (2000). 

[10] G. KhahuUin, Phys. Rev. B 64, 212405 (2001). 

[11] Y. Taguchi, Y. Tokura, T. Arima, and F. Inaba, Phys. Rev. B 48, 511 (1993). 

[12] Y. Okimoto, T. Katsufuji, Y. Okada, T. Arima, and Y. Tokura, Phys. Rev. B 51, 9581 
(1995). 

[13] D.A. McLean, K. Seto and J.E. Greedan, J. Sol State Chem. 40, 241 (1981). 
[14] C. Ulrich, M. Reehuis, unpubhshed data. 

[15] V.G. Zubkov, I.F. Berger, A.M. Artamonova, and G.V. Bazuev, Sov. Phys. CrystaUogr. 
29, 296 (1984). 

[16] E.F. Bertaut, Acta Cryst. A 24, 217 (1968). 

[17] As our data agree very weU with an analysis based on spin canting alone, we have not 

9 



considered effects arising from a possible alternation of the magnetic form factor due to 
orbital order, as discussed in Ref. 



[18] M. Tsubota et al, Physica B 281 & 282, 622 (2000). 

[19] J. Kulda, Program RESTRAX (private communication). 

[20] G.S. Rushbrooke et al. in Phase Transitions and Critical Phenomena, edited by C. Domb 
and M.S. Green (Academic, New York, 1974). 

[21] See, e.g., S.W. Lovesey, J. Phys. C 5, 2769 (1972). 

[22] M. Itoh, M. Tsichiya, H. Tanaka, and K. Motoya, J. Phys. Soc. Jpn. 68, 2783 (1999). 

[23] J. Akimitsu et al, J. Phys. Soc. Jpn. 70, 3475 (2001). 

[24] S. Ishihara, T. Hatakeyama, and S. Maekawa, Phys. Rev. B 65, 064442 (2002). 

[25] G. Khaliullin and S. Okamoto, Phys. Rev. Lett., following paper. 



10 



FIGURES 



0.8 



0.6 - 



^ 0.4 
c 
g 

^- 0.2 - 




zi- 



0.0 - 



c 
CD 

E 0.10 
o 



0.08 - 



g 0.06 
CO 

E 0.04 



0.02 
0.00 



YTiO, 



a) 



(0, 2, 0)o 
(1>1-0)c 




H \ ^ \ ^ 1 " I T I 




(0, 1, 1)o 
(1/2, 1/2, 1/2)c 



10 20 30 40 

Temperature (K) 



50 



FIG. 1. (a) Ferromagnetic and (b) G-type antiferromagnetic components of the ordered mag- 
netic moment of YTiOs extracted from the amphtudes of the (0, 2, 0)o and (0, 1, l)o magnetic Bragg 
reflections, respectively. A weakly temperature dependent nuclear contribution to both reflections 
has been subtracted. The inset gives a pictorial representation of the magnetic structure. The 
TiOe octahedra are highlighted. 
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FIG. 2. Representative constant-q scans with (a) q = (|, j,0)c and (b) q = ^, at 5 K 
(< Tq) and 50 K (> Tq). At 5 K both phonons and sharp magnons are present. The magnons 
become overdamped above Tc- The hnes are results of a convolution of a magnon cross section 
with the spectrometer resolution function, as discussed in the text. 
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FIG. 3. Spin wave dispersion relations at T = 5 K. The closed symbols refer to magnetic 
excitations, whereas the open symbols can be identified as phonons. The grey points were extracted 
from scans at low temperatures, without a corresponding scan above Tc, so that an unambiguous 
identification is not yet possible. The solid line is the magnon dispersion relations derived from an 
isotropic Heisenberg model. The dashed line is the same model augmented by a spin anisotropy 
as discussed in the text. 
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FIG. 4. Exchange constants (in units of Jse) as a function of the Hund coupling r) for different 
values of the parameter n^y characterizing the orbital state. The inset shows the ratio of in-plane 
and out-of-plane exchange couplings for fixed 77 as a function of Uxy 
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